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We study the polymeric nature of quantum matter fields using tiie example of a Friedmann- 
Lemaitre-Robertson- Walker universe sourced by a minimally coupled massless scalar field. The 
model is treated in the symmetry reduced regime via deparametrization techniques, with the scale 
factor playing the role of time. Subsequently the remaining dynamic degrees of freedom are polymer 
quantized. The analysis of the resulting dynamic shows that the big bang singularity is resolved, 
although with the form of the resolution difTering significantly from that of the models with matter 
clocks: dynamically, the singularity is made passable rather than avoided. Furthermore, the results 
of the genuine quantum analysis expose crucial limitations to the so-called effective dynamics in 
loop quantum cosmology when applied outside of the simplest isotropic settings. 



I. INTRODUCTION 

Einstein's theory of general relativity (GR) successfully 
describes gravitational phenomena, predicting with high 
precision all large scale observations made to date. It is 
however expected to fail in the ultraviolet regime due to 
the quantum nature of the reality at the Planck energy 
scales. To obtain accurate predictions for such situations 
one has to resort to quantum gravity (QG). 

Despite many attempts [H-Q , no general, complete and 
working (quantitatively) formulation of QG exists. In 
particular, in the context of the canonical quantization 
programs, so far it was possible to complete the quantiza- 
tion program only in certain situations, where gravity is 
coupled to specific matter fields (for irrotational dust see 
0, 0] and also the earlier, well defined implicit construc- 
tions in order to generalize these frameworks 
(or complete the alternative approaches) it is crucial to 
first study in detail the simplified mini- and midisuper- 
space settings. 

In this paper we consider the simplest minisuperspacc 
model, which represents a Friedmann-Lemaitre-Robcrt- 
son- Walker (FLRW) universe - an isotropic, flat space- 
time admitting a massless scalar field as source. This 
model is widely used as a testing ground for QG methods 
and is at the same time of particular interest in cosmol- 
ogy 

In the context of QG, the model has been studied in 
detail using tools of loop quantum cosmology (LOG - see 
the references in the second paragraph and |lll4l7l) . In 
LQC, an application of the polymer quantization |l8l - [2^ 
to the geometric degrees of freedom results in a dynam- 
ical singularity resolution [l^l, whereby the big bang is 
replaced by a big bounce. This result was later confirmed 
(at the genuine quantum level) for different matter fields, 
in particular the Maxwell field [l^] and dust [13] . 



However, the abovcmcntioncd big bounce result was 
obtained via a somewhat "hybrid" approach: the ge- 
ometry is quantized via loop techniques, while the mat- 
ter (the scalar field) is treated by methods of standard 
quantum mechanics (Schrodinger representation) Q Fur- 
thermore, the system was analyzed by methods dedi- 
cated to theories with a time reparamctrization freedom. 
Namely, the evolution was implicitly defined by means 
of the formalism of partial observables |2^. The above 
approach was also applied outside of the isotropic set- 
tings, both for homogeneous models (like various Bianchi 
models [l^]) and inhomogeneous spacetimes (in particu- 
lar Gowdy models [l^) as well as in the context of per- 
turbation theory about the cosmological sectors [29l - l3]| n 
The results for the inhomogeneous settings are however 
based on heuristic methods and the dynamic is not sys- 
tematically investigated. 

An alternative approach is presented in [32| (see also 
jssj). where the investigated model is the one consid- 
ered here but including a nonvanishing cosmological con- 
stant. The analysis is carried out in the context of quan- 
tum geometrodynamics, which is based on the standard 
Schrodinger quantization of the metric Hamiltonian for- 
mulation of Arnowitt, Dcser, and Misncr (ADM). More 
specifically, the system is treated via the so-called de- 
parametrization technique (for an example starting from 
the full theory see @): one of the dynamical variables 
- in this case the scale factor - is selected as clock at 
the classical level, after which the system can be quan- 
tized and regarded as freely evolving with respect to this 
clock. The obtained results are (by the majority) consis- 
tent with those of the studies of the same systems in the 
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Recently, the results of [17| were confirmed I2a| through the anal- 
ysis of the same model, where both the geometric and the matter 
degrees of freedom are quantized via polymer techniques. This 
work used one of several possible in this context loop quantiza- 
tion schemes (see [iol and the discussion in Sec. Illll l. 
The list of references given here contains only selected examples 
representing the current state of development for each model. 
For a more complete list we refer the reader to [Tsl . ITst . 
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framework of gcomctrodynamics in |34|- l36l | , w here the 
partial variable formalism has been applied [26| . 

A consistent treatment requires the quantization of the 
geometry and the matter in the same way. In the context 
of LQC the intermediate step towards this goal is the 
analysis of a loop quantized scalar field coupled to gravity 
quantized via standard techniques. This step is necessary 
to identify the physical effects arising specifically due to 
the polymer nature of matter. 

In full loop quantum gravity (LQG) a consistent quan- 
tization of the scalar field was proposed in [37| , where one 
of two possible (and inequivalent) implementations of the 
polymer representation [l^ was used. The same choice 
was later made in [2^ to derive the symmetry reduced 
description and to determine the LQC dynamics of (this 
form of) the polymer scalar field. 

The alternative (in a certain sense dual to the above) 
consistent quantization prescription was applied in [s^, 
where again the FLRW isotropic universe is investigated 
by scalar quantum mechanics on a classical cosmologi- 
cal background. Elements of a semiclassical analysis led 
to the construction of an effective approximation of the 
dynamics, of which the study showed that the big bang 
singularity is replaced by a past-eternal de Sitter phase 
("eternal inflation") with graceful exit. 

The mathematical formalism characteristic to this pre- 
scription was later successfully extended to the inhomo- 
geneous setting in the context of quantum field theory 
on Minkowski space and on a cosmological background 
[39} . Both of these extensions were built via Bojowald's 
lattice refinement techniques jiot . 

Since the nature of the studies mentioned just now is 
semi-heuristic, a comparison with the genuine quantum 
dynamic is indispensable. This is exactly the goal of the 
work presented here. To provide the precise quantum 
theory we first perform a deparametrization analogous to 
the one in [32|, choosing the scale factor cubed as clock. 
Then, we quantize the scalar field via loop techniques, 
applying the prescription originally provided in (ssj . 

In our work wc focus on the precise construction of the 
quantum model, that is in particular the correct defini- 
tion of the Hilbert space, and the analysis of the physical 
consequences: the dynamic, the existence of a semiclas- 
sical sector, and the correct GR limit. Surprisingly, the 
requirement of the latter will have a critical impact on 
the form of the Hilbert space and, consequently, on the 
domain of applicability of the heuristic construction of 
the so-called effective dynamics from LQC. 

The numerical analysis of the dynamic shows that, yet 
again, there is no quantum big bang. However, instead 
of bouncing back, the quantum state (the wave packet) 
transits deterministically through the point marking the 
singularity in GR. The quantum evolution picture ap- 
pearing here resembles thus the one advertised in the 
"early LQC epoch" Consequently, this work (see 



also the results in [25|) suggests that the big bounce is 
an effect arising solely due to the polymeric (discrete) 
quantum nature of the geometry. 

The paper is structured as follows: in Sec.|ll]wc in- 
troduce the details of the classical FLRW model that we 
analyze. Then we proceed in Sec. IIIII with the construc- 
tion of the precise quantum theory. Finally, in Sec. lIVI 
we analyze the physical results and conclude in Sec.lVl 
with a general discussion. 

In our studies we select the natural units c = h = 1 
and introduce the abbreviation = 127rGN for a length 
scale. Later in the paper we further restrict our attention 
to the case L = 1 corresponding to a form of Planck units. 



II. CLASSICAL THEORY 

In this paper we focus on the case of an isotropic and 
flat FLRW universe with a minimally coupled massless 
scalar flcld cj) = 4>{T) e M as source. The metric tensor of 
such a universe can be expressed as 



-(iVdr)^ + a^JydXMX^' 



(2.1) 



where N = N{T) e is the lapse function and a = a{T) e 
M is the scale factor. Since the "symmetric criticality 
principle" [1, |4l| is valid in the present situation, the 
canonical action 



A- 



(Pa 

Ti 



a + Pd, 



H[N])dT, 



(2.2) 



can be directly and conveniently derived from the re- 
duced Einstein-Hilbert action. The canonical momenta 
appearing in (j2.2p are 



Pa 

Vc 



/ 3 \^ \a\a 



Vc 



N 

{0,^4 = 1, 



(2.3a) 
(2.3b) 



{a,Pa} = l, 

whereas the scalar Hamiltonian constraint takes the form 

, 2 



H[N] = N 



2X2 



/ L^Pg Y ^ I LP^ Y 



(2.4) 



Note that to arrive at the form (|2.2p of the action we had 
to first introduce the 3 + 1 splitting R x J\f, where for 
the flat universe (considered here) Af = R^. Due to the 
noncompactness of the spatial slices we had to introduce 
in the process of deriving (|2.2p the infrared regulator - 
a cube or "cell" V c finite size (see [H, IH] and 

the discussion in [1^). Its physical volume is V = Vc|ap, 
where Vc = /y d^X is the comoving coordinate volume of 
V. 



^ We note that the early results were derived for a different system, 



where the geometry instead of the matter was polymeric. 
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This step introduces an additional complication into 
the treatment as one has to make sure that the result- 
ing model has a well defined (unambiguous) regulator 
removal limit. The classical FLRW theory is invariant 
under the rescaling X i-^ (X, C e R^, which increases V 
by a factor of C'^ (an active diffcomorphism) . This in- 
variance is a natural requirement for the description of 
the model to remain well defined when the regulator is 
removed [i^. We stress that this requirement is how- 
ever by no means sufficient in the quantum theory (see 
in particular psj). Furthermore, in the quantum the- 
ory the C-invariance is not given trivially [ij, [13, EH and 
hence imposing it as a condition for consistency affects 
the choice of the canonical variables for the quantization 
[see (|2.1ip , p.l2p , and the paragraph prior to them] . 

The first-class |45l - l47j Hamiltonian constraint H[N] « 
generates infinitesimal transformations of T. As can 
be seen from (|2.ip . there also is the possibility of a 
reparametrization of the scale factor a and the Euclidean 
metric 5 by an 77 € K''= such that a >->■ a/rj and S 1-^ 77^ (5, re- 
spectively. This residual 77-symmetry corresponds to the 
freedom of fixing the coordinate scale (a passive difFeo- 
morphism) and orientation (a "large gauge transforma- 
tion" [i3,|l3,|33,li2,lil). Just like it is required for the 
(■-transformation, the physics has to be invariant under 
an r^-transformation. Among the ry-invariant quantities 
are the action A, the volume V and ratios of the scale 
factor such as the Hubble parameter h = a/{Na). 

In the next step we fix the time reparametrization free- 
dom by implementing the seconrf-class [45[447| gauge 



G = T- 



Vca^ 
L2 



(2.5) 



We emphasize that the time-gauge G becomes T - 
sig{r])Vca^ / under an 77-transformation. This means 
that the orientation of a relative to T changes if ry e . 
Replacing T by -T has no impact on the space of so- 
lutions to the Wheeler-DeWitt equation resulting from 
(12.71) . and thus amounts to a "time-reversal" operation 
[49|. However, the reduced classical and quantum for- 
malism based on (j2.5|l and derived here is the result of 
singling out one of ±r, and is therefore no longer time- 
reversal invariant. If (|2.5p defines a future-directed clock, 
the past orientation would be given by the gauge con- 
straint G = T + Vca'^/L^. These considerations are rele- 
vant in the construction of the initial state for the quan- 
tum evolution (see Subsec. llVBp . 

The canonical formalism we constructed here has the 
deficiency of explicitly depending on the infra-red regula- 
tor V since, according to (|2.5p . the clock variable T scales 
like C"^ . This would make the removal of the infra-red reg- 
ulator V from the resulting quantum theory a rather te- 
dious task. As the initial step in addressing this problem 
we replace T by the dimensionless variable t = T/\T'^\, 
where T'^ € is some fixed but otherwise arbitrary ref- 
erence value. Furthermore, since also scales like 
we analogously define = P4,/\P^\ with P^ e being a 
fixed reference value for P^ . Technically, the replacement 
of T by f can be brought about by a simple change of the 
integration variable in (j2.8p . whereas the replacement of 
Pip by Pcj, is realized by an "extended canonical" or "scale 
transformation" [H^l- Altogether, this procedure yields 
the canonical action 



A='-^ {Lp^cf>'-Hs) 

L JlT^lti 

where (/>' = d(f)/{dt), {(f), pc/,} = 1/L, and 



Given the equation of motion 



dt, 



(2.11) 



{a,H[N]} = -N 



9Vc\a\ 



(2.6) 



the form of G implies in particular that Pa is negative. 
This and the form of the constraint H[N] then lead to 
the reduced canonical action 



A= r^^(p^0-i/R,)dT, 

JTi 



(2.7) 



where again {cf^P^} = 1. The reduced Hamiltonian takes 
the form 



-Pt 



L'Pa 



L'Pn 



3Vca2 



3Vca2 

Finally, the consistency condition 

dTG + {G,H[N]} = {) 
uniquely determines the lapse function 



P^ 



LT 



(2.8) 



(2.9) 



(2.10) 









t 


1 


t 



dt- 



■Hr dT, 



(2.12) 



is the Hamiltonian related to Hji by a scaling. We stress 
that t and arc dimensionless variables. 

At this point it is necessary to mention that the above 
modification does not yet remove the dependence of the 
theory on the infrared regulator. Indeed, while the con- 
stant T"^ and P^ are fixed, no particular value of V can be 
distinguished on physical grounds. In consequence, the 
particular "physical" universe is represented by classes of 
solutions rather than by single ones. This nonuniqueness 
can be easily shown on the level of specifying the initial 
data. There, the single universe regulated by different 
cells V will correspond to the entire set (class of equiv- 
alence) of the initial data at the chosen initial time t* . 
This dependence will propagate through to the quantum 
theory. 

For initial t* ,p*^ e the solutions to Hamilton's equa- 
tions of motion derived from (|2.1ip and (|2.12p are 



m 



(2.13) 
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Therefore, the lapse function A'' = \T^\l\P'^Piji\ is a con- 
stant [positive since we fixed the time orientation by the 
chosen gauge constraint p.Sp ], The canonical represen- 
tation of the spacctime Ricci scalar takes then the form 
R = -6h^, which in turn implies 



3R{t) 



(^) ■ = ■ <"^) 



The form of the reference value i?'' suggests the nat- 
ural and simplifying choice \P^\ = |T''|, which corre- 
sponds to fixing the initial curvature value to be R{t*) = 
-{2/3)[p'^/{Lt*)]'^ . This partially fixes the freedom of 
rescaling t and but still does not remove the implicit 
^-dependence discussed earlier. 

In order to simplify the expressions, from now on we 
will use the Planck units normalized by i = 1. With this 
choice the energy density g and the pressure p are 



Q{t)=p{t) 



-Im, 



(2.15) 



and the spacetime singularity occurs at t = 0*. The pos- 
itivity of N implies that the future-pointing evolution of 
the scalar field is "into" a big crunch for t e ST and "away 
from" a big bang for t € M"^ . From (|2.13|) it is evident that 
for an element of the branch of the solution space admit- 
ting a big crunch, the value of (j){t) - </>* for e is 
related to the analogous value for € by an overall 
sign-change. The same holds for for an element of the 
big bang branch of the solution space. Furthermore, we 
have the large time-inversion symmetry 

(sigipimt) - = -(sig(p;)[0(t) - 0*]),, (2.16) 

which relates big crunches (the left-hand side for < e R") 
with big bangs (the right-hand side for i e R^). In the 
canonical formalism at hand these identities are manifes- 
tations of the invariance of the covariant action under a 
replacement of (j) with -(f>. 

Finally, we note again that the variables t and p^ 
are dimensionless but still not invariant under a 
transformation. The observable scalar field in (|2.13|) and 
the spacctime Ricci scalar in (|2.14p - along with its re- 
lated scalars in ()2.15p - are inheriting this implicit non- 
invariance. This fact will play a crucial role in singling 
out the correct regularization scheme in the quantum the- 
ory. 



III. QUANTUM THEORY 

Our goal here is to build the precise quantum mechan- 
ical representation of the model introduced above. This 
means in particular the construction of a suitable Hilbert 
space H and the representation of Hg as a self-adjoint op- 
erator acting on the suitable domain in H . The quantum 
evolution will then be determined by some Schrodinger 



equation 



for ip = '4'{t, I 



I — w 



i' = Hsip 



(3.1) 



A. The scalar field momentum operator 

To start with, let us recall that the canonical formal- 
ism introduced in the previous section describes a freely 
evolving, isotropic and fiat FLRW model. The evolution 
is governed by the Hamiltonian Hg in (|2.12p . which by 
((231) depends on the scale factor clock t = Vba^/|T''|. The 
scalar field (j> is thus the only object subject to a quanti- 
zation. Here, we have several possibilities to proceed. 

The most obvious way to construct the quantum de- 
scription is to apply the Schrodinger representation, as 
in j32l |. As shown there, this representation does not 
lead to a singularity resolution, as the semiclassical wave 
packets simply follow the classical trajectories into the 
singularity. 

An alternative approach, which is pursued here, is the 
implementation of the polymer representation. The re- 
quirement of the existence of an infra-red regulator re- 
moval limit (see the previous section) forces this repre- 
sentation to be time- dependent^ 

To begin the detailed specification of the polymer 
quantization procedure, let us briefly recall the stan- 
dard Schrodinger representation. It is characterize d by 
the Stone- von Neumann uniqueness theorem p^.[20l[5^. 
which implies that among all the irreducible regular re- 
alizations of the Wcyl form 



hJu 



of the canonical commutation relation 



[<i>,P4>] = ii 



(3.2) 



(3.3) 



on the space L^(R, d0) of Lebesgue square-integrable 
functions, the Schrodinger representation 



h = e' 



Jfj, = e 



(3.4) 



is unique. The regularity property says that the map- 
pings of A to I\ and /x to are continuous, which holds 
if for tl',uj 6 L^(M,d(/)) the mappings 



are continuous. 



(3.5) 



* This situation is analogous to the one in the loop quantization 
of the geometric degrees of freedom, where consistency require- 
ments label the improved dynamics construction as the correct 
one [H. 
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To generalize the above formalism, let us now con- 
sider the space of exponentiated operator labels, further 
parametrized in the following way 



(3.6) 



The subscript t can be seen as parametrizing the param- 
eters A and fi of the groups of unitary operators Ix and 
Jp, respectively. That is, the Weyl algebra given in p.2p 
depends now on time and so do the unitary operators de- 
fined in p.4p . However, we observe that because of the 
continuity the operators (f) and p^ are 



= i lini — -^-^ = i lirn — -^-^ = 01, 



At^O At 



>y^o At 



lim — = -i lim ■ 



(3.7a) 



ilTT, (3.7b) 



thus they are independent of t. To conclude, in the 
Schrodinger quantization the operators defined in p.7p 
are not changing if the group parameters are themselves 
parametrized by t. 

The situation changes drastically if the regularity con- 
dition in the Stone- von Neumann uniqueness theorem is 
dropped. In this case a possible faithful realization of 
p.2|) is given by the so-called "polymer representation" . 
For the sake of generality we will further allow it to be 
time- dependent (in a yet unspecified way). The non- 
separable Hilbert space H for this representation con- 
sists of functions ip e P{M.,^^) satisfying the square- 
summation requirement 



The inner product on H providing this norm is 



(3.8) 



(3.9) 



where uj is another element of P(M, ^^). We denote by 
^0 the measure that maps a subset of M to its cardinality 
(the so-called "counting measure") and by 

Dm^W^UL^^C-^o), B^,(V,^) = ULm*('/'o), (3.10) 

(/.oeS(i/')/s 0o6(S(V')nS(i.-j))/s 

domains defined in terms of the necessarily countable 
support S of '0 or w. For 0, % 6 M we have the equiv- 
alence (/) S X if and only if there is an integer i e Z such 
that (j)= x + ipt [recall p.6p ] . The domains D^^ are then 
disjoint unions of uniform lattices 

h^AM = {M + ^Pt, 0oe[O,M*)- (3.11) 
Orthonormal basis states of H are "half-deltas" 



1, '/' = X, 

0, otherwise. 



(3.12) 



extending the definition of the Kroneckcr delta symbol 
to the real line. 

The time-dependent polymer representation is now 
given by 



(3.13) 



which characterizes again a multiplication and a trans- 
lation operator, respectively. The "A-mapping" given in 
p.5[) is once more continuous for At so that by Stone's 
theorem [l^ [s^] the scalar field multiplication opera- 
tor remains to be given by p.7p . The difference to 
the Schrodinger representation in p.4p is that the "/z- 
mapping" in p.Sp is no longer continuous in fif There 
is therefore no self-adjoint momentum operator gener- 
ating infinitesimal translations. On P{M.,^^) there is 
only an operator generating finite translations. We are 
therefore forced to regularize it, for which we employ the 
technique introduced by Thiemann in the context of full 
LQG [s^l. In essence this technique is approximating 
the undefined Pri, by well-defined translation operators. 
Following [13, [33, [sj] we choose 



^Mt 




(3.14a) 
(3.14b) 



The action of the former on a state ?/; e H is 



[i^{(f> + Pt)-H<l^-Pt)] (3.15) 



so that, if we could send pt to [or according to 
p.6p send p to 0, thereby taking the limit at the kine- 
matic level], we would get back the differential opera- 
tor -\d/{d(j))- We observe that the representation of the 
momentum operator p^p,, is highly non-unique, in the 
same way the representation of finite difference opera- 
tors in numerical analysis is. We stress that, unlike in 
the Schrodinger representation, because of p.6p the mo- 
mentum operator is now time-dependent. 

At this point it is necessary to emphasize that the pre- 
sented polymer quantization is not the only possible one. 
Essentially, by replacing the roles of (j) and we arrive 
at another polymer representation, inequivalent (and in 
a sense "dual" ) to ours (see the discussion in [l^ ) . Such 
a dual representation was used in the quantization of 
the scalar field in full LQG [13]. Its symmetry reduced 
version was applied to the LQC model of an FLRW uni- 
verse [2^ filled with a massless scalar field. The subse- 
quent analysis of the spectral decomposition of the evolu- 
tion operator (playing the role of the Hamiltonian) shows 
that the dynamic of such a system is exactly the same as 
the one of the system with the scalar field quantized via 
standard methods of quantum mechanics |34l ] . Both ap- 
proaches, ours and the one of [25j . are equally viable from 
a mathematical point of view. Therefore, choosing one 
of them requires a physical input. 
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B. The Hamiltonian 

The next step is the construction (including the deter- 
mination of the action) of the quantum Hamiltonian Hs 
generating the unitary evolution. To do so, we switch 
to the scalar field momentum space, which is again the 
Pontryagin dual of the real line but this time the latter 
is equipped with the discrete topology. In short, it is 
the Bohr-compactified real line Kb- The Hilbert space 
defined in the previous subsection is then equivalent to 
the space H\ which consists of Bohr square- measurable 
functions 

^\Pc^) = EV'(0)e-'*''* 6 L2(RB,(dp^)B) (3.16) 
satisfying 



1 r'^' 



B 



^^^^ Ij^'^p^^^' 



< oo. 



(3.17) 



The inner product (between and another w*! € H*!) 
generating this norm is 

{i;^\u;^) ^ f i^oj^p^) {dp^)B. (3.18) 

The basis orthonormal with respect to it is formed by the 
plane waves 



(3.19) 



We observe that for a uniform lattice I}f^^{ip) = L^j(0o) 
[see p. lip ] the Bohr measure {dp^)B becomes the 
Lebesgue measure dp^ with an integration over the fixed 
interval (-7r//Xt, 7r//.t(]. The momentum space polymer 
theory defined here would then be that of Fourier with 
discreteness in position rather than momentum space. 

In the general polymer theory at hand, the action of 
the multiplication and translation operator on the basis 
states 60 is unchanged in comparison to (|3.13p so that 
the scalar field operator is given by \d/{dp^). On the 
other hand, the operator is now the regularized p^^j 
(see the previous subsection), which is approximated by 
translation operators. Indeed, since e^+^j = e^^e^ we 
have 



Jfit ~ ^ P<}>fj.t 



sin(^fP0) ^ 



(3.20) 



because of which the action of the Hamiltonian operator 
can be explicitly given by 



sin{fitp^) 



p-tt 



(3.21) 



where i/j^ = ip^{t,p^). Given this, the fact that ||'0|| = \\ip^\\ 
allows us to immediately prove the conservation of the 



norm under an action of Hs- To show this we write 
explicitly the time derivative of the norm 

5„,,„o r.d 



= E7^(i,0)V(t,x) f iJs6-'f*(*->^) (dp0)l 



(3.22) 



To evaluate the right-hand side we first observe that the 
integral can be expressed as the contour integral 



lim 

c^oo 2C 



L r = hn 

:C J 

-c 

i-l)-Klfj.t 

f * f 



-(2i-l)ir/^t 



lim V I f 

™ 4n7r fii \ J 



(2j-l)7r//it 



2iTTl fit 

2iTT/fJ.t 

f 

-(2i-l)7r//it 2{i-l)Tr/tit (2i-1)tt / tj.t 



(3.23) 



The specific form of this integral allows us the to drop 
the absolute value Hs in p.2ip ) replacing it instead with 
a sign appropriate for each integration domain. Next, 
wc apply some simple trigonometric identities, the ^t- 
translation invariance of . 



])^j(-!/'), and (see |55|) 
' 2mr , 



^ cos 

4=1 



■(2i + l)7r 



!sin|— -((^-X)! 



(3.24) 



Finally, if we divide this by n and take the limit n ->• oo 
[sec diini)], wc get i9||?/'||/at = 0. 

Up to now, the shift function /xj in the approximated 
p^ operator is an arbitrary function of time. At the math- 
ematical level the situation is analogous to the one in the 
loop quantization of the geometry (see jSGj). where the 
fiducial holonomy length could be an arbitrary function 
on the phase space). There, however, the physical consis- 
tency requirements restricted the possible choices to just 
one class of functions . We expect that the same situ- 
ation occurs in our model. To show that this expectation 
is indeed realized let us recall the following facts. 

The particular moment of the universe's evolution can 
be represented by various points on the phase space cor- 
responding to different choices of the regulator cell. Fur- 
thermore, once we ask about the locally measurable prop- 
erties of the universe (observables) at this moment, there 
has to exist their nontrivial limit as we remove the regu- 
lator. 

One such local observable is the energy density (|2.15p 
determined by (j2.14p . The quantum operator corre- 
sponding to it is related to the Hamiltonian Hs in the 
following way 



1 



-Hi 



(3.25) 



From (|3.2ip it follows that at the fixed time t the spec- 
trum of this operator equals 



Sp(e) 



0, 



1 



2M?i2 



(3.26) 
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The most natural way to satisfy the consistency require- 
ments discussed in the previous paragraph is to require 
that the spectrum of g be time independent. This imphes 
/it oc 1/t so that we can fix the hmction pt in (j3.6p as 



Pt 



(3.27) 



This in turn gives Ht = p/\t\, which for the "volume clock" 
V = t/p results in the momentum space Schrodinger equa- 
tion 



ov 



(3.28) 



for = ■0'i(u,p^). 

Note, that the above method of fixing fit is almost a 
full analog of the conditions used for the geometry de- 
grees of freedom in [sij . There, however, the reasoning 
exploited the existence of "nicely" behaving semiclassical 
sectors through the use of the so called effective dynam- 
ics. Here, as we have not yet investigated the dynami- 
cal sector implementing that reasoning directly would be 
risky. Instead, we managed to fix the ambiguity through 
the considerations on the genuine quantum level. 

In the next section we solve the Schrodinger equation 
p.28p in order to analyze the dynamics of the system and 
to discuss the physical properties of the solutions. 



IV. THE DYNAMIC 

The Hilbert space and the explicit action of the Hamil- 
tonian operator constructed just now allow us to easily 
determine the system's dynamic. At this level the re- 
quirement of the theory to be physically meaningful be- 
comes crucial. The principal requirement is that the the- 
ory must have the proper low energy limit. Here, this 
means that in the distant past and future the quantum 
evolution ought to agree with the predictions of GR. In 
our case an inability of the model-description to realize 
this property would imply that the formulation should 
be further and adequately corrected. In fact, as we will 
see below, this is precisely what is required here. 

To begin, let us investigate the dynamic of the theory 
exactly as specified in the previous section. 



A. Single Lattice Hilbert space 

Once we select as the configuration variable, the 
Schrodinger equation given in (I3.28P becomes a simple 
ordinary differential equation, which is easy to solve on 



the domains t > and t < 0. Its solution reads 

i>\v,P4,) = Eyy*1p\v* ,P4,) 



S{v,p^) 



sm 



p<f> 




V 


) V _ 



(4.1a) 

(4.1b) 
(4.1c) 



where we set v* 
gral function is 



t* /p. For |arg(2;)| < tt the cosine intc- 
cos{y) 



Ci(z) = 7 + ln(z) + 



Jo 



1 



y 



dy (4.2) 



with 7 the Euler-Mascheroni number [53| . The definition 
implies Ci{z) 7 + ln(z) for z ^ 0, suggesting semiclas- 
sical behavior of sufficiently sharply peaked initial states 
in the limit 00. 

At the point 1; = the operator Hs is not well defined 
so that the Cauchy-Peano theorem, which normally en- 
sures the existence and uniqueness of the solution, cannot 
be applied there. Thus, in principle one cannot extend 
the solution to (|3.28p through that point. However, the 
solutions (|4.ip admit well defined limits v 0"^ for v > 
and V < 0, and thus one can define an extension to that 
point by taking these limits. Since 



lim F(v, 



= 0, 



(4.3) 



there exists a unique unitary operator evolving states to 
the instant v = 0. Furthermore, it has the very simple 
form 



(4.4) 



As a consequence there exists a preferred extension of 
the evolution through w = 0, defined by the requirement 
of continuity of ■0'' at w = 0. The global solution is again 
given by (|4.ip . 

It appears that the existence of such a preferred ex- 
tension is sufficient for singularity resolution. However, 
as we will see below, this is not the case. To explain 
what is missing, we consider any normalized initial state 
il^{v* ,p^) such that the expectation value of the scalar 
field operator is finite 



< 00. 



(4.5) 



Then the expectation value of (j) at any volume v is given 
by the following formula 



cl){v) = {i;\v\4>\'iP\v) 



-S{v*,P4,\ 



5(^;,p,)Ci(|^|) 



Since the cosine integral function belongs to L^( 
the above expectation value is in fact equal to (/ 



(4.6) 

^,dp0), 
. That 
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is to say the evolution is frozen. This result is then in 
direct disagreement with the predictions of GR. In con- 
sequence, our states exhibit an unphysical behavior in 
the low energy (large v) limit. 

Our model then still lacks an appropriate physical 
Hilbert space. To explore the possibilities of construct- 
ing it, let us first go back to analyzing the solutions to 
p.28p but this time by considering the wave functions on 
the configuration space as opposed to the momentum one 
used in (|4.ip . Then, the evolution of a state tpy = tpiv,-) 
can be viewed as an assignment w i-^ e H^,, where w e M. 
The Hilbert space H„ is spanned by eigenstates of Hs for 
a fixed value of v. However, per analogy with the loop 
quantization of the geometry [3J] we can distinguish sec- 
tors that are invariant with respect to the action of Hg at 
V. These sectors consist of functions that are supported 
on the lattices L((/?o)a*«j where fiy = l/\v\ and 

U^o)=^^AM/l^v = Wo} + ^, <^oe[0,l). (4.7) 

We can then attempt to consider at the initial v = v* the 
subspaces ^vcpo - ^v\fj.^L{^g) as the superselection sec- 
tors and evolve them independently. Such a decomposi- 
tion can be performed at each v independently. Let us 
now probe whether there exists any relation between the 
spaces H^j^Q for different values of v. The answer is given 
immediately by the form of the solution (|4.ip : since the 
cosine integral function is non-periodic, the unitary evo- 
lution to any v + v* immediately couples an infinite num- 
ber of these lattices. In consequence, the sectors H^;^,-, 
of hL, are not true superselection sectors in the sense of 
[H [ll m, m, 13 . Therefore, we are forced to work with 
the original non-separable Hilbert space H„ without ac- 
cess to previously available tools that allow to distinguish 
separable subspaces. Furthermore, the form of (|4.6p sug- 
gests that in order to provide a nontrivial evolution, the 
physical Hilbert space needs to be equipped with a con- 
tinuous rather than a discrete inner product. 

B. Integral Hilbert space 

A similar situation appeared in LQC already in a differ- 
ent context during the studies of the FLRW universe with 
a massless scalar field and a positive cosmological con- 
stant [1^. There, following the choice of a lapse adopted 
to using the scalar field as time variable, the evolution op- 
erator admitted a family of self-adjoint extensions, each 
with a discrete spectrum. However, a different choice 
of the lapse - corresponding to parametrizing the evolu- 
tion by the cosmic time variable - led to a unique self- 
ad.ioint generator of the evolution with continuous spec- 
trum jSq . The physical Hilbert space corresponding to 
the latter case appeared, furthermore, to be an integral of 
all the Hilbert spaces corresponding to the particular self- 
adjoint extensions of the former case, with the Lebesgue 
measure determined by the group averaging procedure. 

Motivated by this observation, we introduce the ana- 
log of the integral Hilbert space in our case. First, we 



note that on the domain [0, 1) of ipo one can introduce 
a natural (quite general and time dependent) Lebesgue 
measure M{v,ipQ)dipo. Next, we introduce a decompo- 
sition of the non-separable Hilbert space H onto spaces 
H^ipo at the initial time v* . This defines the decomposi- 
tion of the initial data at v = v* 

H 3 V(f * , 0) 1^ ip^o =Tp{v\ (j>) , (</,o) e ^v*vo ■ 

(4.8) 

This initial data is then extended to the solutions to 
(|3.28p via (|4.ip . We thus have a decomposition of the 
physical Hilbert space onto explicitly separable (at least 
at V = V*) subspaces. 

Now, we compose the new Hilbert space via an integral 

Hp„. = £ H,*^„M(v*,(^o) d(^o, (4.9) 
equipping it with the inner product 

(4.10) 

where ipvipo = V'¥'o(^i<^) ^ ^vipo- This is our candidate for 
the physical inner product: between each pair of solutions 
it is evaluated on the initial data slice at v = v* . On 
that initial slice it can be written in a very simple form, 
namely, 



(V'«*|a;„*) = J^tp{v*,(j))u}{v* ,(t))M{v* ,(poi(t)))/Hv* dcj). 

,,(4.11) 

It is by definition time-independent but a priori it may 
not have a simple local form analogous to (|4.1ip atv + v*, 
which can potentially complicate the evaluations of the 
expectation values of the observables. 

We note, however, that the construction performed for 
V = V* can be repeated at each value of w, giving rise to 
potentially inequivalent constructions of the candidate 
Hilbert space. One can then consider a function 

P{ipy\ujy) = / 'ilj{v,(t))uj{v,(t))M{v,ipo{<f)))/i-Lvd(l). 
Jr 

(4.12) 

On each slice of constant v this function equals the inner 
product of the candidate Hilbert space constructed with 
respect to this slice. One can then ask under which con- 
dition these Hilbert spaces will be equivalent and their 
inner products equal. A condition necessary and suffi- 
cient for it is that dP{ipy\uj.u)/dv = 0. The form of the 
unitary evolution operator (|4.1ap implies however that 
this condition will be satisfied if and only if 

M{v,ipo) = ^j.ym{ipo). (4.13) 

Following this choice, our candidate Hilbert space be- 
comes (up to a rescaling ip^ >->■ ifiyl^/m on H„) the space 
L^(R, d0) with the standard L^-inner product. Also, the 
momentum space is now L^(R, dp^) with the correspond- 
ing Lebesgue measure. 
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Using (|4.4p . we can now consider an initial state 

^\v\p^) = E,.^i^\Q,p^) (4.14) 
with a real.^ L^-normalized Gaussian 




-i"^(p«-pJ)^/2 



(4.15) 



This class of states is "special" in the sense that the 
quantum evolution they undergo is semiclassical both for 
V ~ V* and v — v* (see below). Furthermore, the states 
-E'^oV'HO?^^) with unit L^-normalized ■(/'''(Ojp^) also span 
the solution space of the Wheeler-DeWitt equation de- 
fined by the Hamiltonian constraint in (|2.4p and with 
Vca^ = T identified as clock. They are thus particularly 
convenient in comparing the evolution in the Schrodinger 
and polymer quantizations. The set of the complex con- 
jugate of these states represents the analogous states for 
the clock -Vco?' = T (see [i^ and recall that we set 
L=l). 



We can now try to evaluate once more the expectation 
value of the scalar field. For that we chose an example of 
a "semiclassical" Gaussian state peaked at v* = 250 with 
w = 1. This choice results in a relatively small initial 
value of both the scalar field and its momentum fiuc- 
tuations (see below). We further set = 5 to prevent 
any significant portion of the Gaussian initial state from 
overlapping with the momentum space origin, as this is 
where the cosine integral function has an integrable sin- 
gularity. The quantum evolution of such initial state was 
then calculated numerically. Fig. [1] shows the quantum 
trajectory corresponding to this evolution. From there, 
it is evident that the evolution is semiclassical for |f | » 1. 
In fact, the classical solution that is well approximating 
the quantum trajectory is characterized by 



+ H\vl) 



(4.16) 



where the overall sign "=f" corresponds to v e R=^ and, 
where v* is taken to be the largest plotted value of v (250 
in the present case) . Given the definition of the function 
S, we observe that |?;*| » Ip^I/tt is a necessary require- 
ment for semiclassicality. As we can see, the physical 
state is indeed passing in a continuous manner through 
the point v = 0, corresponding in the classical theory to 
the big bang singularity [which is particularly clear from 
Fig. 1(a) and also from (14. 3p ]. This happens regardless 



of the sign of the initial momentum so that the quan- 
tum evolution is effectively respecting (|2.16p . which in 
the classical theory specifies the relation between the so- 
lutions for negative and positive t. 



-H 

St 
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V = t/p 

(a) Quantum evolution for two values of . 



Ill 













-5 



-10 -5 5 10 
V = t/p 

(b) Classical and quantum evolution for = 5. 

FIG. 1. An illustration of the quantum evolution of the scalar 
field is presented in Fig. 



1(a) In Fig. 1(b) we display the 



= 5 quantum evolution along with the corresponding clas- 



sical evolution over a subset of the i;-interval used in Fig. 1(a) 
The circles, triangles and squares correspond to actual mea- 
surements made using the software Matlab and Octave. 
The time-interval between two consecutive measurements be- 
comes smaller as \v\ approaches 0. In this plot, the smallest 
value of \v\ is 0.25 but values as small as 10"'* have been consid- 
ered with the same outcome. Namely, there is semiclassicality 
for \v\ » 1 and the big bang singularity at ?; = is resolved. 



To examine more closely the issue of the singularity res- 
olution we also analyzed the expectation values of the op- 
erator corresponding to the spacetime Ricci scalar. The 
quantum trajectory is presented in Fig.O The measure- 
ments are independent of the overall sign of p^, thus only 
the case p*^> has been plotted. We see that the space- 
time curvature remains finite for all v. This confirms the 
analytical result of (|3.26p and, thus, implies the global 
boundedness of the spectrum of the Ricci scalar operator 
once pv is fixed via (|3.27p . One can thus conclude that 
the big bang singularity is resolved. 
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<^ -1 




-200 





V = t/p 



200 



FIG. 2. Expectation values of the spacetime Ricci operator 
for various values of v are plotted in this figure. As in Fig.[T] 
measurements are taken the more often the closer \v\ gets to 
and the smallest value is |t;| = 0.25. It follows that the big 
bang curvature singularity is resolved and the branches ii e 
are connected. 



0.720 



0.715 - 



0.710 



-200 





v = t/p 



200 



FIG. 3. In this figure an illustration of the fluctuations of the 
scalar field operator is presented. For \v\ -» the fluctuations 
increase but remain nonetheless finite. 



Finally, in Figs. [3] and |4] the fluctuations of the scalar 
field and polymer momentum operator are depicted. 
They both quickly approach a constant value (for the 
presented example approximately equal to 2"^^^ for both 
quantum fluctuations) as \v\ increases, which is the value 
expected for a Gaussian with w = 1. This confirms the 
semiclassical nature of the state for \v\ ~ v* . What is 
interesting in the near-singularity region is the fact that 
the fluctuations of the scalar field operator are in fact 
decreasing for 0. This may happen because the 

state gets "squeezed" towards the origin in order to "fit 
through" the point v = 0. This however requires a more 
detailed analysis of the nature of the state there, which 
may be the subject of a subsequent investigation. 



ii- 

^ 2 




-200 





V = t/p 



200 



FIG. 4. This figure provides an illustration of the fluctuations 
of the polymeric momentum operator. Just like for the scalar 
field fluctuations in Fig.|3] the behavior for \v\ » 1 indicates 
semiclassicality. 



V. DISCUSSION 

In this paper we investigated the quantum dynamics 
of the isotropic and fiat FLRW universe of infinite ex- 
tent and sourccd by a minimally coupled massless scalar 
field. Our focus was on the modifications to the dynam- 
ics following from the polymeric nature of the matter and 
in particular the issue of the singularity resolution. To 
single out these effects we implemented one of two pos- 
sible loop quantization schemes of the scalar field. This 
scheme is the analog of the one used so far in LQC to 
quantize the geometry degrees of freedom. Unlike in the 
most previous works in LQC, instead of implementing the 
Dirac program to solve the Hamiltonian constraint, we 
performed a complete deparametrization of the system 
by choosing a scale factor dependent time variable. As a 
result, the physical evolution is described by a free Hamil- 
tonian. The quantization of such a deparametrized sys- 
tem is implicitly equivalent to selecting the Schrodinger 
quantization for the geometry when applying the Dirac 
program. Therefore, the effects of the geometry discrete- 
ness are not featured in our model. 

In the process of constructing the correct description 
of the quantum system we encountered several obstacles: 

First, the noncompactness of the universe's spatial 
slices forced us to introduce an infrared regulator. The 
necessary consistency condition that the theory has to 
admit a well defined and nontrivial regulator removal 
limit restricted then the Hamiltonian to particular form, 
which happened to be explicitly time-dependent. 

Second, the Hilbcrt space to which the physical states 
belong occurred to be non-separable. This is a standard 
(and treatable) problem in LQC. Here however, the ex- 
plicit time dependence of the Hamiltonian prevented us 
from implementing the known technique of subdividing 
the (too big) Hilbcrt space onto separable superselection 
sectors. An idea to naively proceed with determining 
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the dynamics on that space led to a model significantly 
disagreeing with GR predictions at the low curvature 
limit. Indeed the quantum evolution of the scalar field 
was frozen. 

To cure this defect we performed a specific construction 
of the separable Hilbert space out of the nonseparable 
one, taking as the guideline the relation between Hilbert 
spaces corresponding to the models with different choices 
of the lapse function in LQC in the presence of a posi- 
tive cosmological constant. As a result, we were able to 
construct a certain integral Hilbert space equipped with 
continuous rather than discrete (as usually in LQC) inner 
product. 

Such a construction of the Hilbert space was then used 
to investigate the dynamics. To do so we selected a class 
of Gaussian initial states and evolved them numerically. 
The resulting quantum trajectories showed a good con- 
vergence to the classical trajectories predicted by GR at 
low energies. At high curvatures (small \v\) however we 
observed a significant departure from GR. Indeed, the 
most critical feature of the model is the existence of a 
unique unitary evolution operator evolving to/from the 
time slice v = corresponding to the classical singularity. 
This and the regularity of the wave function describing 
the physical state allowed us to select a naturally pre- 
ferred extension of the evolution, thus ensuring a deter- 
ministic evolution through the classical singularity. Fur- 
thermore, the quantum counterparts of the Ricci scalar, 
energy density, or pressure are explicitly bounded oper- 
ators. In consequence, the listed quantities remain finite 
throughout the entire evolution, including in particular 
v = 0. 

At this moment it is important to note that, unlike 
in previous contributions to the literature on this model, 
here the quantum features responsible for singularity res- 
olution originate from the matter rather than the geomet- 
ric sector. Therefore, the form of the singularity resolu- 



tion is different than in the literature: Instead of being 
avoided, the surface w = is made passable and all the 
standard locally measurable quantities remain finite. 

Finally, let us comment on an important lesson learned 
from this model: the predicted dynamics depends criti- 
cally on the construction of the physical Hilbert space 
of the model, even though the regularized form of the 
Hamiltonian remains the same. This implies in particu- 
lar that the regularized form of the classical Hamiltonian 
or Hamiltonian constraint is not sufficient to robustly 
determine or even well approximate the quantum evolu- 
tion. This issue is particularly critical in all the studies 
of models in LQC performed via the so called effective 
dynamics techniques without prior specification of the 
elements of the genuine quantum system. 

In a further project we intend to take a closer look 
at the behavior of the state near the singularity. Why 
do the quantum fluctuations of the scalar field decrease 
towards the origin of the time-axis? Of interest is also 
the inclusion of a non-zero cosmological constant. Fi- 
nally, and this is most intriguing, we would like to address 
the question of how the quantization procedure presented 
here can be combined with that of the geometric sector 
discussed in the LQC works [H, [13, [M El El ■ 
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